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Abstract: Following Campos's work [Phys. Rev. Lett. 88, 141602 (2002)], we investigate 
the effects of temperature on flat, de Sitter (dS), and anti-de Sitter (AdS) thick branes in 
five-dimensional (5D) warped spacetime, and on the fermion (quasi-)localization. First, in 
the case of flat brane, when the critical temperature reaches, the solution of the background 
scalar field and the warp factor is not unique. So the thickness of the flat thick brane is 
uncertain at the critical value of the temperature parameter, which is found to be lower 
than the one in flat 5D spacetime. The mass spectra of the fermion Kaluza-Klein (KK) 
modes are continuous, and there is a series of fermion resonances. The number and lifetime 
of the resonances are finite and increase with the temperature parameter, but the mass of 
the resonances decreases with the temperature parameter. Second, in the case of dS brane, 
we do not find such a critical value of the temperature parameter. The mass spectra of 
the fermion KK modes are also continuous, and there is a series of fermion resonances. 
The effects of temperature on resonance number, lifetime, and mass are the same with the 
case of flat brane. Last, in the case of AdS brane, the critical value of the temperature 
parameter can less or greater than the one in the flat 5D spacetime. The spectra of fermion 
KK modes are discrete, and the mass of fermion KK modes does not decrease monotonically 
with increasing temperature parameter. 
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1 Introduction 

In the braneworld scenario our universe is a 3-brane embedded in a higher dimensional 
spacetime. The braneworld scenario has received a lot of attentions since it can provide 
us a novel approach to resolve the cosmological constant and the hierarchy problems [1-4], 
and reproduce the Newtonian law of gravity [5]. In the braneworld scenario all kinds of 
matter fields should be localized on the brane, and there is much literature on these issues 
[6-32]. The braneworld scenario also provides some new approaches to resolve the family 
and favor problems in particle physics [33-39]. 

Further, taking into account that our real universe is a system with the temperature. 
So it is interesting to investigate the effects of temperature in braneworld theories. Brevik et 
al. [40] investigated the quantum (in)stability of the AdSs braneworld universe at nonzero 
temperature. Campos discussed the critical phenomena of flat thick branes in warped 
spacetimes [41] with a five-dimensional (5D) complex background scalar field. Bazeia et 
al. [42] investigated the geometric transitions of thick braneworlds at high temperature 
limit. Different from the ideal thin braneworld models [4, 5], thick braneworld models 
usually need to introduce 5D background scalar fields to generate the branes [41, 43-69]. 
To introduce the temperature effects to thick braneworld theories, we can calculate the 
effective potential of the 5D background scalar fields at finite temperature. And the method 
is the same as the one proposed in [70-72] in four-dimensional fiat spacetime. Ansari and 
Suresh [73] calculated such an effective potential of cj)'^ model in 5D fiat spacetime with the 
one-loop correction. As far as we know, the effective potential of 5D scalar fields in thick 
braneworld theories has not been calculated. To do such a calculation is very meaningful 
but also difficult, not only because the spacetime is curved, but also because the background 
scalar field is coupled with the gravity. Although without an analysis function form of the 
effective potential with the temperature, one can still investigate the temperature effects 
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in braneworld theories [41, 42]. Within the (j)'^ model, only the mass parameter is corrected 
by the effects of temperature, so the variation of the mass parameter will reflect the effects 
of temperature on branes. 

Our study is based on the work of Campos [41]. In this paper, the self-interaction 
potential of the background complex scalar field ^ has the form of [41] 

= a|$|2 - bM^R - 30?) + + Cs, (1.1) 

where and 0i are the real and imaginary parts of the complex scalar field respectively, 
and C5 is a constant. If there is no background complex scalar field C5 will be the 
cosmological constant of 5D spacetime. Taking the effects of temperature into account, the 
mass parameter a should be a function of the temperature T [70-72]. But the parameters 
b and c will not vary with the temperature. This can be seen from V{^) after a shift of 
$ around its vacuum expectation value ^o(0ROi </'io): 'I' — )• f^^I' + [(j)K — >• + (/)roi 
—7- 6cl)i + (pio]. After that shift the coefficient before |$o|^ is also c, so c will not vary 
with the temperature. The coefficients before (j)'^^ and <?i>Ro0io changed from —b and 36 
to —b + 2C50R and 36 + 4c5(/)r, respectively. Because the vacuum expectation value of 64>: 
(O|50r|O) = 0, so the coefficients b will also not vary with the temperature. Because we 
have not the exact function form of a(T), we only take a as the temperature parameter in 
this paper. 

Campos found that the presence of gravity would lower the critical value of the tem- 
perature parameter of the phase transition comparing to 5D flat spacetime [41] . And the 
phase transition is characterized by the emergence of a double kink solution of the (pi. 
Below the critical value of the temperature parameter, the solution of the (pi has a single 
kink form. Further, we find that, in the case of AdS brane, the critical value of a can less or 
greater than the one in the 5D flat spacetime, and we do not find the double kink solution 
of the (pi for the case of dS brane. 

We also investigate the effects of temperature on fermion localization and quasi- 
localization on flat, dS, and AdS thick Branes. In order to localize fermions on thick 
branes, the coupling between fermions and the background scalar fields is needed. We 
consider the Yukawa coupling i]^(pi^. We find that, in the cases of fiat and dS branes, 
the fermion zero mode can be localized on the branes, and there are quasi-localized mas- 
sive fermion KK modes, namely, fermion resonances. The number and the lifetime of the 
fermion resonances increase with the temperature parameter a. But the resonance mass 
decreases with the increase of a. In the case of AdS brane, the spectrum of fermion KK 
modes is discrete, and the number of the discrete KK modes is infinite. The variation of 
masses of fermion KK modes does not decrease monotonically with the increase of a. 

The paper is organized as follows. Firstly, the effects of temperature on flat, dS, 
and AdS thick branes are discussed in section 2. Secondly, the effects of temperature 
on the fermion localization and fermion resonances are investigated in section 3. For 
self-completeness and self-consistency, the Schrodinger equations and the corresponding 
potentials for the KK modes of left- and right-handed fermions are derived. Before per- 
forming numerical calculations, we analyze the asymptotic behaviors of the potentials at 
zero and infinity along the extra dimension. In sbusection 3.1, we investigate the effects 
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of temperature on the number, lifetime, and mass of fermion resonances quasi-localized on 
flat and dS branes. We show the variation of masses of fermion KK modes localized on 
AdS brane with the temperature parameter a in subsection 3.2. Finally, the conclusion 
and discussion are given in section 4. 

2 Effects of temperature on thick branes 

To investigate the effects of temperature on thick branes, we consider a 5D action with a 
complex scalar field minimally coupled to gravity: 

S = j d'xdy^ Qi? - ^g''''dM<f*dN^ - V{^)^ , (2.1) 

where V{^) is the same as the one in Eq. (1.1). To ensure the stability, the parameter c 
must be positive. There are three degenerate global minima for the potential V{^) if 6 > 
and < a < 9ac/8, where ac = iP'/^c [41]. The three degenerate global minima are at 
= ,^0, = {-\ + i4)^o, and q>3 = (-^- where <Ao = 1(1 + Jl - ^) [41]. 



The line element in this model is assumed as 

= e^^^'\g^,^{x)dx^'dx'' + dz^), (2.2) 

where iijf = 0,1,2,3, ^^jy is the 4-dimensional metric on branes, and e^"^^^-* is the warp 
factor. The usual hypothesis is that A and $ are only the functions of the extra dimension 
coordinate z. 

For the three types of maximally 4-symmetric branes (flat, dS and AdS), the equations 
of motion of A{z), (pi, and have a unified form: 

^'l = -3A'cP[ + e'^^, (2.3) 

4 = -3^% + e2^^, (2.4) 



^" = ^"-^((0f + 0g) + A4), (2.5) 

6A'' = 2A4 + ^((0f + <P'i) - 2e'^V), (2.6) 

where the prime stands for the derivative with respect to z. A4 is the cosmological constant 
on branes, and A4 = 0, > and < correspond to the cases of flat, dS and AdS branes, 
respectively. Equations (2.3)-(2.6) are a set of nonlinear ordinary differential equations, 
which can be solved numerically ^ with the following boundary conditions: 

A{0) = A'{0) = (t^'niO) = 4>m = 0, (2.7) 

1 a/3 

(/)r(+oo) = --00, and (pi{+oo) = —(j)o. (2.8) 



^We can solve these equations with the FORTRAN code colsys and bvp_solver, or the function bvp5c in 
Matlab. 
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Figure 1. The profiles of the background scalar fields 0i and (/>r for flat brane. The three lines 
correspond to a = 0.5, 0.8, 0.961, respectively. The other parameters are set to 6 = 2 and c = 1. 




Figure 2. The profiles of the warp factor A for flat brane. The three lines correspond to a = 
0.5,0.8,0.961, respectively. The other parameters are set to & = 2 and c = 1. 



The solutions of A, and are varied with the temperature parameter a. For flat, 
dS, and AdS thick branes the variation of the solutions of (pi and (/>r with a is similar. 
So, we only show the solutions for the case of flat brane in Fig. 1 with 6 = 2, c = 1, and 
different values of a. Fig. 1 shows that the solutions of flat thick brane are dilated along 
the extra dimension with the increase of a, and the solution of cpi tends to have a double 
kink form. 

The solutions of A{z) for flat and dS branes are similar, and the solutions for the case 
of flat brane are shown in Fig. 2. 

For the case of AdS brane, the behavior of the warp factor is very different from the 
cases of flat and dS branes. This has been analyzed in Refs. [29, 74, 75]. In this case, with 
the metric (2.2) and the potential (1.1), the z coordinate only runs from —zt, to Zj,, where 
Zfe is a finite value, and A{z) is divergent at itz^. This can be seen from the analysis of the 
solutions of Eqs. (2.3)-(2.6) at z — )• iz^. When z — )• izj,, 

1 a/3 

M±Zb) = --cPo, M±Zb) = ±^<Po, (2.9) 
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Figure 3. The profiles of tlie warp factor A{z) for AdS brane. The three lines correspond to 
a — 1.07, 0.9, and 0.5. The other parameters are set to 6 = 2, c = 1, and A4 — —0.1. 



and 



b'^i±Zb) = cP[i±Zb) ^ 0. 



So Eqs. (2.5) and (2.6) turn into 



A"{z ^ ±Zb) = A'\z ^ ±zh) - -A4, 
6A'\z ^ ±Zh) = 2A4 - e^^^^^^^^V,,, 



(2.10) 

(2.11) 
(2.12) 



where Vz^=V{4>i{±Zh), (j)ji{±Zh)) is a constant. The solution of the above equations is 



A{z ±Zb) log , 



'2A4 



log 



-A, 



sm 



(2.13) 



So it is clear that the warp factor A(z) is divergent at the boundaries z = ±Zb- This is 
consistent with the numerical solution of A{z) shown in Fig. 3, from which it can be seen 
that the value of Zb decreases with increasing a. 

The critical value of the bulk temperature parameter a is defined such that the profile 
of (pi has the form of a double kink [41], namely 4>i{0) = 0. For the case of five-dimensional 
flat spacetime {A = 0), we can easily determine the critical value of a, which is given by 
Oc = 6^/(4c) [41]. But for five-dimensional warped spacetime with flat brane, the critical 
value of a is not Oc but a smaller effective critical value a* [41]. 

We find when a = o^, the solutions of cpi, (p^, and A are not unique ^. This means that 
at the critical temperature the width of flat thick brane is not fixed. The profiles of (pi and 
(pn with a* = 0.9617 (b = 2, c = 1) are shown in Fig. 4 for flat brane. In the AdS brane 
case, we find that the effective critical value a* can greater than ac- But for the dS brane, 
we do not find a double kink solution of (pi. 

In the Appendix we have discussed the relation between a* and Oc with the assumptions 
that as a —7- a*, |(ac — a*)/ac| <C 1 and (/)r(0)^ <^ (pQ. And we show that, for the case of 



^ Because Eqs. (2.3-2.6) are non-linear, so in order to solve them numerically we need to ofTer initial 
guess solutions of 0i, ^r, and A. And different guess solutions follow different solutions of them. 
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Figure 4. The profiles of (jii and with different solutions with a■^, = 0.9617 for flat branes. The 
other parameters are set to 6 = 2 and c = 1. 




Figure 5. The profiles of (f>i and (/)-[{ with a = 1.08, a ~ 1.09, and a = 1.10 for AdS brane. The 
critical value a* is 1.084. The other parameters are set to 6 = 2, c = 1, and A4 = —0.1. 

flat brane, a* < Oc- But in the cases of dS and AdS branes, the relation between a* and Oc 
is uncertain. 

In the flat brane case, when the value of a is greater than the critical value a*, we do 
not flnd any solution of Eqs. (2.3)-(2.6). However, in AdS brane case, we do, and we show 
the proflles of (pi in Fig. 5 for o = 1.08, a = 1.09, and a = 1.10 (the critical value a* is 
1.084). 

3 Effects of temperature on fermion localization and resonances 

In order to localize a bulk fermion on these branes, we need to introduce the coupling of 
the fermion and the background scalar $. A general form of the coupling is r]^F{(f)i, 
where ^ is the fermion field and F{(j)i, (j)^) is a function of and i^r. The action of the 
fermion describing such coupling is 

S^= f d''xV^[^r''{dM+uJM)'^ - r?^'F(0i,0R)^}, (3.1) 
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where = (e ^7^,e ^7^) are the curved space gamma matrices, rj is the couphng 
constant and 77 > 0. ojm is the spin connection [76-78] and its nonvanishing components 
are [14, 79] 

where cD^ is the spin connection derived from the metric c/fj^u- The resulting Dirac equation 
is 

7^(a^ + u^) + 7'(5. + 2A') - r,e^F{(t>i, </.r)] ^ = 0, (3.2) 

where 7^(9^ + '^At) is the Dirac operator on branes. With the usual KK and chiral decom- 
position 

^{X, Z) = Yl (V'Ln(x)/Ln(^) + V'Rn(x)/R„(z)) , (3.3) 

n 

where subscripts 'L' and 'R' denote the Left- and right-handed chiralities, respectively, and 
the four-dimensional massive Dirac equation 

7''(9^ + a;^)^Ln = "T-n^Rn, (3.4) 
'~i^{d^+Qj^)'ll)'R,n = mnlp-Ln, (3.5) 

we can obtain the coupled equations of /Ln and /r„ [23]: 

[dz + rie^F{(t)i,(t)K)\fhn{z) = m„/R„(z), (3.6) 
[dz - ??e^i^(0i,</'R)]/Rn(^;) = -mnfhniz), (3.7) 

where /Ln and /r„ satisfy the following orthonormality conditions: 

/oo 
fLmjRndz = (^mn^LR- (3.8) 
-oo 

Further, we obtain the Schrodinger equations [23]: 

[ - + C/fL(z)] /l„ = mlhn, (3.9a) 
[ - 52 + C/fR(z)] /r„ = mlfj,n, (3.9b) 

where the effective potentials are given by 

UiLiz) = V^e^^F{cl)i, 0r)2 - r]e^F{cPi, ./.r)' - 7]F{<Pi, 0R)e^A', (3.10a) 
UfKiz) = UfLiz)\r,-,-^. (3.10b) 

In order to determine the forms of Ufi,{z) and C/fR(z), we should first construct the 
form of F((f)i,<f)ji). There are two simplest forms of F: Fi = (pi and F2 = (/>r. With the 
solutions of (pi and (pn, we find that (pi is an odd function, and (pn is an even function. So, 
in order to hold the Z2 symmetry of Uji^{z) and U^{z), we choose the first form of F, i.e., 
Fi = (Pi. So 

C/fL(z) = 7fe^^(pf - r]e^(p'i - v^ie^M, (3.11a) 
Uik{z) = v'^e^^(Pf + r?e^0'i + r](Pie^A'. (3.11b) 
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Figure 6. The profiles of the potentials of fermion KK modes for flat brane. The parameters are 
set to 77 = 2, 6 = 2, c = 1, and a = 0.5, 0.8, 0.961. 

3.1 Flat and dS branes 

Before solving Eq. (3.9) numerically, let us first analyze the behavior of U^{z) and U{r{z) 
at z = and z — t- ±00. 

First, at 2; = 0, with the boundary conditions A{0) = A'{0) = (/)i(0) = 0, we get 

UiL{0) = (3.12a) 
Um{0) = (3.12b) 

And from Fig. 1, we see that (/>j(0) > 0, so [/fL(0) > and f/fR(0) < for positive coupling 
constant 77. 

Second, when z — >■ ±00, (j)[{z) = 0. So we have 

UfLiz ±00) = r]'^e^'^<l)f - 7?(/)ie^^', (3.13a) 
Um.{z ±00) = ry^e^^c/^i^ + 7?(/.ie^^'. (3.13b) 

And Eqs. (2.5) and (2.6) turn into 

A" = _ Ia^^ (3.14) 
3 

6A'^ = 2A4 - e^^Foo, (3.15) 

where Voo=V {(j)i{±oo) , 4)ii{±oo)) is a constant, and (/>r(±oo) = — |(/>o, (/)i(±oo) = it-^c^o- 
In the flat brane case, A4 = 0, the solution of A{z — )• ±00) is 

A{z ^ ±00) = - log (V-Ko/6|z| + ci) , (3.16) 

where ci is an integration constant. So e'^'^\z^±oo = l/(\/— Voo/6|2;| + ci)^\z^±oo — ^ 
and A'{z)\z^±oo = -l/i^/-Voo/6\z\ + ci)| 

z—^±oo ~^ 0- Then we reach the conclusion that 

f^fR — ^ and UfL — ^ when z — )• ±00. 
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Figure 7. The profiles of the zero mode of left-handed fermions for flat brane. The parameters 
are set to ?] — 2, b — 2, c — 1, and a = 0.5, 0.8, 0.961. 



In the dS brane case, A4 > 0, the solutions of A{z), e^^^\ and A'(z) at z 
A'iz) = ^/f tanh (^^i\z\ + c,)^ 




±00 are 
(3.17) 

(3.18) 

(3.19) 



Then we come to the conclusion that and f/fR vanish at 2; — ?• ±00. 

With the numerical solutions of A{z), (pi, and (pn, we plot the potentials C/fL and C/fR 
in Fig. 6 for flat brane. The profiles of the potentials for dS brane are similar with the case 
of flat brane, and we do not show them again. We can see that, for rj > 0, the potential 
of left-handed KK modes f/fL is a modified volcano-type potential, which has a well lower 
than zero, so only the left-handed fermions have zero mode. The Schrodinger Eqs. (3.9) 
can be solved using the Numerov method [80] with two initial conditions at z = 0, which 
can be set as 



/L,R(0)=cil,4R(0)=0 

for the even parity KK modes and 

/L,R(0)=0,//R(0) = d2 



(3.20) 



(3.21) 



for the odd parity KK modes, where di and d2 are arbitrary constants. Here we choose 
di = d2 = 1- The profiles of the zero mode for flat brane are shown in Fig. 7 with different 
values of the temperature parameter. The case for dS brane are similar. 

The volcano-type potential implies that there will not exist discrete spectrum of 
fermion KK modes, but there may exist fermion resonances [79, 81-83]. The massive 
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Figure 8. The profiles of Pl.r for left-handed fermions for flat brane with a = 0.961, 6 = 2, c = 1, 
and 1] = 2. 



KK modes cannot be normalized because their wave functions are oscillating when far 
away from the brane along the extra dimension. So Ref. [82] proposed a function 



I. 



/_\l/L,R(^)Prfz 



(3.22) 



as the relative probability for finding the resonances on a brane. Here 2zc is about the width 
of the thick brane and z^ax is set to z^^ax — Wzc- So for KK modes with > I/^X 
{Vj^^ is the maximum value of Vl,r), /l,r can be approximately taken as plane waves, 
and the value of Ph,R{rn) will trend to 0.1. 

As an example, we plot the profiles of PL,R{m') in Fig. 8 corresponding to a = 0.961 
for flat brane. In this figure each peak corresponds to a resonant state. We can estimate 
the lifetime r of a resonance as r ~ where F = dm is the full width at half maximum 
of the peak [84]. It can be seen that the spectra of resonances are the same for both the 
left- and right-handed fermions. 

In order to investigate the effects of temperature on resonances. We solve Eq. (3.9) 
with the following values of a: 

a = {0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 0.92, 0.94, 0.96, 0.961, 0.9612, 0.9614, 0.9616}. 



-2 < Vf:''^ into account. 



Here we only take those resonances with in 

Because the spectra of resonances or bound states for left- and right-handed fermions 
are the same, here and after we only discuss the spectrum of left-handed fermions. 

In order to see the effect of temperature on the fermion potentials and Ufi^ more 
intuitively, we can take the volcano box potential approximation [45] for the potentials. 
For example in Fig. 9, the well's depth, height, and width are Vi, V2, and 2zi, respectively. 
The width of barrier is \z2 — zi\. Now back to see the potential C/fL in Fig. 6, we can find 
that with the increase of a the width of barrier and the width of the well are increasing, but 
the depth and height of the well are decreasing. Because the variation of height of the well 
is small and all the KK modes with > 0, from the knowledge of quantum mechanics, we 
can expect that with the increase of a the number of the resonances is increasing and their 
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Figure 9. The fermion potential C/fL with a = and its volcano box potential approximation. 
The other parameters are set to 6 = 2, and c—\ 
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Figure 10. The variation of the resonance mass m with a for left-handed fcrmions for flat brane. 
The other parameters are set to = 2, c = 1, and rj ^ 2. 

masses are decreasing. The numerical results support this expectation. The variation of 
the resonance mass m with a is shown in Fig. 10. We can see the resonance mass decreases 
with the increase of a. While the resonance lifetime r increases with a, which is plotted in 
Fig. 11. 

The variation of the number of resonances n with the temperature parameter a is 
shown in Fig. 12. It can be seen that the number of resonances n is increasing with the 
temperature parameter a, and there is a threshold value for a to produce more resonances. 
For example, for a < 0.4 there is no resonance, and for 0.4 < a < 0.9 there is a resonance. 
This is very similar to the phenomenon of the particle generation in colliders, namely, if we 
want to find more particles we need to increase the center-mass energy and the generation 
of particles need the center-mass energy greater than a threshold value. 

3.2 AdS brane 

When z — )■ 0, the behaviors of the potentials of left- and right-handed fermions for AdS 
brane are similar to the ones in the flat braneworld scenario, so we will not repeat them 
here. 
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Figure 11. The variation of the resonanee hfetime r with a for left-handed fermions for flat brane. 
The other parameters are set to = 2, c = 1, and rj = 2. 
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Figure 12. The variation of the number of resonances n with a for left-handed fermions for flat 
brane. The other parameters are set to 6 = 2, c = 1, and rj ~ 2. 

But when z — t- ±Zf,, the behaviors of A are very different from the ones in the cases 
of flat and dS branes, so we need to analyze the behaviors of the potentials of left- and 
right-handed fermions carefully. When z — t- izf,, the potentials turn into 

C/fL = rie^Mve^(^i - A'), (3.23a) 
UfR = rje^U7]e^^i + A'). (3.23b) 



Inserting Eq. (2.13) into the above equations, we get 

-A4 77 {rj/rjo - cos {z')) 
3 rjo sin^(z') 



^fL = -i^z "72777^ ' (3.24a) 



= Z^:L (^/^o+;of)) ^ 

3 r/o sin^(z') 

where z' = y^-A^/Sizb - \z\), r]o = with vq = V{(j)i = ^(/)o,<Ar = -^M- Then we 
can see that, if 77 > ijo, 

U[l{z ±zb) +00, (3.25a) 

Um.{z ±2b) +00; (3.25b) 
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Figure 13. The plots of U[l and J7fR for AdS brane. The other parameters are set to = 2, 6 = 2, 
c = 1, and A4 = —0.1. 



ifr] = r]o, 



UiLiz^±Zb)^-^, (3.26a) 
U{r{z ±Zb) +00; (3.26b) 



if < r/ < 



C/fL(z — )• ±Zb) — )• —00, (3.27a) 
C/fR(z ±Zb) +00. (3.27b) 

In the cases oi t] = 7]q and < < r^o, the behaviors of left- and right-hand fermions 
potentials at z — )• zizZb are so different, and we will not expect the left- and right-handed 
fermions have the same KK modes spectra. So we only discuss the case of r] > rjQ. The 
plots of the potentials of left- and right-hand fermions are shown in Fig. 13. 

From Fig. 13 and the above analysis, we know that, for rj > ijq, all massive fermion 
KK modes are bound states and there is a discrete mass spectrum for them. 

We use the method proposed in Ref. [85] to calculate the fermion mass spectrum. And 
we calculate the fermion mass spectra with a set of values of a: 

a = {0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1, 1.08, 1.10}, 

and the other parameters are set to = 2, 6 = 2, c = 1, and A4 = —0.1. The masses of 39 
KK modes are obtained. 

From Figure 13, we can see that the profile of potentials ?7fL and Uf^i is similar to that 

of the infinite square-well potential. And we know that a particle in the infinite square- 

2 

well potential will hold the the energy of En oc^ [n = l,2,3,...], where d is the width of 
well, and En decreases with d. Figure 13 shows that with the increase of a the width of 
the potentials is decreasing, so we except that the masses of the KK modes increase with 
a. But this is not all consistent to our numerical results. For the 28th excited state and 
the ones above it, the masses of KK modes increase with a. For example, we show the 
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Figure 14. The to — a curves of 28th-35th excited states for left-handed fermfon with parameters 
77 = 2, & = 2, c = 1, and A4 = -0.1. 
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Figure 15. The m — a curves of lst-12th excited states for left-handed fermions for AdS brane 
with parameters 77 = 2, 6 = 2, c = 1, and A4 = —0.1. 

m — a curves of 28th-35th excited states in Fig. 14. But the ones of below the 28th excited 
state masses of the KK modes do not vary monotonically with the increase of a (see Fig. 
15). This is because with the increase of a that the values of the bottom of potentials are 
also not varying monotonically with the increase of o. For example the first excited state 
(left-handed fermion) mass mi{a = 0.2) greater than mi(a = 1.1) is because the value of 
the bottom of C/fL at a = 0.2 is greater than the one of C/fL at a = 1.1. 

4 Conclusion and discussion 

With the model used in this paper, we found the effects of temperature on flat and dS 
branes are similar, except the existence of the critical temperature. In the case of flat 
brane, there is a critical temperature. While it does not exist for dS brane. The thickness 
of flat brane is uncertain at the critical temperature. We think this is related to the 
potential given in (1.1). If we choose the (j)'^ model of a real scalar field, there will not be 
such a critical temperature even for flat brane. 

In the case of AdS brane, we found that the critical temperature is greater than the 
one in 5D flat spacetime. This is just opposite to the case of flat brane, for which the 
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critical temperature is less than the one in 5D flat spacetime [41]. 

The variation of the temperature parameter will affect the spectrum of fermion KK 
modes. First, in the cases of flat and dS branes, the spectrum is continuous. We found that, 
for some parameters, there are fermion resonances on the brane. The fermion resonances 
can be thought as quasi-localized fermions with finite lifetime staying on the brane. The 
number of fermion resonances is finite and increases with the temperature parameter a, so 
does the resonance lifetime. This means that the brane with a bigger value of temperature 
parameter can trap more quasi- localized fermions on it. The masses of the quasi- localized 
fermions decrease with the increase of a. Second, in the case of AdS brane, if the Yukawa 
coupling is larger than some critical coupling constant, then there is a discrete spectrum of 
fermion KK modes, and all the KK modes (four-dimensional fermions) are localized on the 
brane. But the masses of fermion KK modes do not vary monotonically with the increase 
of the temperature parameter. 

From the model discussed in this paper we can see: (a) In the case of AdS brane, 
the fermions will be bounded on our brane permanently and the number of fermions is 
infinite. But the model with an AdS brane is unreal. This is because we know that in 
our real universe the 4D cosmological constant is positive, (b) Corresponding to the AdS 
brane, the flat and dS branes are more similar to our real universe. In these cases the 
fermions are quasi-localized states, except the zero mode, with a finite lifetime staying on 
the branes. This means that the fermions can dissipate into the fifth dimension just like 
the quasi-localized gauge field [12] and the quasi- localized gravitons [84, 86-88]. From the 
point of view of 4D observer, one will have the chance to find energy non-conservation in 
the collider. 

If we assume that the temperature parameter a is a monotone increasing function of 
temperature T, then, for a collider experiment, the number of the produced new fermion 
resonances will increase with the center-of-mass energy of beams. This is consistent to 
results of the experiment. Further, our another result that the lifetime of fermions increase 
with the temperature, would be examined if the collider experiment find the signal of 
fermions dissipating into the fifth dimension. 
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5 Appendix 



Here we give the discussion about the relation between a* and Oc- 

First at the point of z = 0, substituting ^(0) = A'{0) = cpiiO) = cp'^iO) = into Eq. 
(2.6) we get: 

2A4 + lcj)[iO) -1/0 = 0, (5.1) 

where 

Vo = y((Ai(o),0R(o)) 

= a4(0) - 64(0) + c4(0) + C5. (5.2) 

Because (j)R{z) is an even function of z, we can Taylor expand (f>R{z) around the point of 
z = 0: 

(pKiz) = do + d2Z^ + d^z'^ + ■■■ . (5.3) 
So (^r(O) = do. Substituting this into Eq. (5.2) we get 

Vo = adl-bdl + cdl + C5. (5.4) 



And Eq. (5.1) is rewritten as 



2A4 + ^(AUO) = adl - bdl + cdl + Cg. (5.5) 

Further, when z tends to its right boundary ( for flat and dS branes it is 00 but for 
AdS brane it is z;,, here we use symbol 00 to stand for it). 

Ko = a|$(cx))|2 - 6</.r(oo)(0r(oo)2 - 30i(oo)2) + c|$(oo)|^ + C5. (5.6) 

Substituting Eq. (2.8) into the above equation leads to 

Voo = acfi - h(t)l + cct>t + C5, (5.7) 

where </.o = |(l + ^(9-8a)K/3). So 

C5 = Ko - {a4>l - h4>l + c<Pl). (5.8) 

Substituting the above expression for C5 into Eq. (5.5) gives 

When (/'j(O) — )• 0, a — )• a*. We assume that dp ^ (/)q (this assumption is supported by 

the numerical results which can be seen from Figs. 1, 4, and 5). Neglecting the terms 
containing square and higher powers of do, Eq. (5.9) reads 

a, = {h4>o-ccg)-^{2ki-V^). (5.10) 
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Further we take a* = ac(l — 5) [cc = 6^/(4c) = 1] and assume that \6\ <C 1. So (po = 
1^(1 + The Taylor expansions of Eq. (5.10) respect to 6 is 

a, = ae(l - ^2 + 45^) _ ^(2A4 - yoo)(l - 2<5 + 75^ - 326^) + (5.11) 
After neglecting the terms containing square and higher powers of 5, Eq. (5.11) reads 

a*=ac- -^{2 Ai-Voo) {I -26). (5.12) 

So in the case of flat brane, we have A4 = and Voo < (this can be found from Eq. 
(3.16)), which leads to a^, < ac- 

But in the cases of dS and AdS branes we cannot confirm whether a* > ac or a* < ac- 
That is because that, in the case of dS brane A4 > and Voo > (this can be found from 
Eq. (3.17)), and in the case of AdS brane A4 < and Voo < (this can be found from Eq. 
(2.13)), we are not able to determine whether (2A4 — Voo) > or < 0. Even in the case of 
dS brane, we cannot confirm the existence of the a* using the numerical method. 

In the above discussion we have made use of two assumptions: \6\ <^ 1 and <^ (pQ. 
This will weak the efficiency of our conclusion. So here we make use of the numerical 
method to check the consistency between the two assumptions and Eq. (5.12). From Eq. 
(5.12) we have 

' = '-b^ + Ac4^.-Voo)- ^'-''^ 

Using the numerical method, we have obtained a set of values of V^, do and 6 corresponding 
to a = a* . All results are shown in Table 1 for the case of AdS brane (the corresponding 
solutions of (pi are shown in Fig. 16) and Table 2 for the case of flat brane (the corresponding 
solutions of are shown in Fig. 17), which show that the value of Oc — a* is almost the 
same with the 6 obtained from Eq. (5.12). So the consistency is satisfied. 

Table 1. The numerical results of Voo, and 6 corresponding to a = a* for the case of AdS brane 
with 6 = 2, c = 1, and Ur ~ 1. 



Voo 


do 


a* 




6 


-0.041001 


-0.000327 


0.999 


0.001 


0.000999 


-0.052821 


-0.000350 


1.028 


-0.028 


-0.0279383 


-0.071759 


-0.000366 


1.051 


-0.051 


-0.0506862 


-0.095345 


-0.000339 


1.070 


-0.070 


-0.0691242 


-0.123884 


-0.000382 


1.084 


-0.084 


-0.082387 
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Figure 16. The solutions of 0i with different values of A4 (given in Table 1) corresponding to 
a = a* for the case of AdS brane. The other parameters are set to 6 = 2 and c = 1.. 

Table 2. The numerical results of Voo, do and 5 corresponding to a = a* for the case of flat brane 
with 6 = 2, c = 1, and Uc = 1. The solutions I, II, and III of 0i are showed in Fig. (17), and they 
correspond to the same value of a* = 0.9617. 
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do 
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solution I 
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